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Flat bands as a route to high-temperature 
superconductivity in graphite 


T.T. Heikkila and G.E. Volovik 


Abstract Superconductivity is traditionally viewed as a low-temperature phenomenon. 
Within the BCS theory this is understood to result from the fact that the pairing of 
electrons takes place only close to the usually two-dimensional Fermi surface resid¬ 
ing at a finite chemical potential. Because of this, the critical temperature is expo¬ 
nentially suppressed compared to the microscopic energy scales. On the other hand, 
pairing electrons around a dispersionless (flat) energy band leads to very strong su¬ 
perconductivity, with a mean-field critical temperature linearly proportional to the 
microscopic coupling constant. The prize to be paid is that flat bands can generally 
be generated only on surfaces and interfaces, where high-temperature superconduc¬ 
tivity would show up. The flat-band character and the low dimensionality also mean 
that despite the high critical temperature such a superconducting state would be sub¬ 
ject to strong fluctuations. Here we discuss the topological and non-topological flat 
bands discussed in different systems, and show that graphite is a good candidate for 
showing high-temperature flat-band interface superconductivity. 

The purpose of this chapter is to propose a route to increasing the critical temper¬ 
ature of superconductivity by searching for special electronic dispersion that would 
promote the superconducting strength. We first show that a huge increase in the 
(mean-field) critical temperature is possible, if a dispersionless energy spectrum, a 
flat band can be created in the system in the absence of the interaction leading to su¬ 
perconducting correlations. We then discuss a few known schemes to generate such 
(approximate or exact) flat bands. 
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1 Superconductivity: pairing energy vs. dispersion 


Within the BCS mean-field theory, the occurrence of Cooper pairing at zero temper¬ 
ature can be studied via the free energy density for the pairing energy A: 

where d is the dimensionality, g <0 describes the interaction strength, and Ep{A) = 
+A'^ is the quasiparticle excitation energy at momentum value p, evaluated in 
a system with the normal-state dispersion £p. The first term in Eq. 0 demonstrates 
that the formation of a gap A decreases the energy of quasiparticles, which fill the 
negative energy levels of Dirac vacuum. The second term is the cost of the formation 
of the gap, which perturbs the vacuum. For simplicity, we consider spinless fermions 
and the gap that does not depend on momentum. 

Requiring A to minimize F^, we get the self-consistency relation 


4 = M / _^ (2) 

2 J {2nnY Ej,{A)’ ^ ’ 

or 

1 = M / _^ (3) 

2 J {27thy Ep{A) ■ 

Equation Q dictates the behavior of A at different dimensionalities d and for dif¬ 
ferent normal-state energy spectra Cp. 

For i-wave superconductivity in conventional metals with an isotropic Fermi sur¬ 
face and dispersion e = vp{p — pf) expanded around the Fermi energy Ep = vppp, 
the integral in Eq. Q is concentrated in the vicinity of the Fermi surface 


1 = 


\g\Ad de 

{2nhYvp Jo y'e'^ + A^ 
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(4) 


Here Vp = density of states in the normal metal; is the area of 

the -dimensional Fermi surface; and Euv <C Sp is the ultraviolet cut-off of the log¬ 
arithmically diverging integral, such as the Debye temperature. This leads to the 
exponentially suppressed gap: 


4 


EuvCXp 



(5) 


and correspondingly to the exponential suppression of the transition temperature T^. 

Situation drastically changes when the spectrum of the normal state has a flat 
band - a region in momentum p where Cp = 0. Since within the flat band Ep{A) = A, 
Eq.Q becomes 
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where is the volume of the flat band in momentum space. Hence instead of the 
usual exponentially suppressed behavior in Eq. Q we have the gap A that is linearly 
proportional to the interaction strength: 


4 = 


\8\Vd 

2{2TihY ’ 


(7) 


Since the critical temperature of the superconductors is typically of the same order 
of magnitude as the gap at zero temperature, the resulting superconductivity may 
exist at high temperatures. 

It is instructive to consider the intermediate case when the quasiparticle spectrum 
is 

/ \M 

Then for M > d,Eq. (|^ gives the power law dependence of transition temperature 
on the coupling constant: 

A . (9) 

In the limit of large M:s> 1, the spectrum ([^ transforms to the flat band concentrated 
at p < po, and the gap (j^ asymptotically approaches the linear dependence on the 
coupling g in Eq. Q. The case with d — 1 and M = 2, where A g^, has been 
considered by Kopaev ID and Kopaev-Rusinov El. 

In the sections below, we consider different systems where the exact or approxi¬ 
mate flat bands could be realized. 


2 Flat band induced by interaction 
2.1 Landau phenomenology 

We start with the flat bands induced by interaction between the fermions. As was 
found by Khodel and Shaginyan 13, the interaction may lead to the merging of 
different fermionic energy levels, which results in the formation of a dispersionless 
band, see also Refs. aa. The effect of the merging of discrete energy levels due 
to interaction has been reported in a recent paper. Ref. Q. 

Their argument is based on the phenomenological Landau’s consideration on the 
derivation of the distribution function Wp of fermions at T =0. It is determined by 
the energy functional S’{n{p)}, whose variational derivative gp is the quasiparticle 
energy. The variation of this functional gives the equation for tip and Cp: 

5^ {n(p)} = J ep5np = 0 . 


( 10 ) 
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8£{n(/5)l = Je(p)5n{p)if*p = 0 

two solutions: 8(p) =0 or 6n{p)=0 



e{p) =0 

Fermi surface 



flat band 


Fig. 1 Illustration of the formation of the Khodel-Shaginyan flat band due to interaction between 
fennions. (Top): the Landau model of the Fermi liquid considers the energy (S'{n{p)} as the func¬ 
tional of the distribution function n(p) of quasiparticles. The variation of the functional gives two 
types of solutions: Ep = 0 or 5/ip = 0. (Lefty, distribution of quasiparticles in the class of Fermi 
liquids. Two regions with solutions Srip = 0 are separated by the surface with solution Ep = 0. 
This is the Fermi surface. Topological stability of this class is demonstrated in Fig.[^ (Right): For 
strong interaction (large interaction constant), the intermediate region with the solution £p = 0, 
may become finite. This is the flat band. Topological structure of such flat band is demonstrated in 
Fig.[3] 


Since the quasipaiticle distribution function is constrained by the Pauli principle 
0 < Wp < 1, there are two classes of solutions of the variational problem. One class 
is gp = 0, which is valid if 0 < rip < 1; another one is 5np = 0 with iip = 0 or iip = 1. 
In the following we illustrate these two classes of solutions. 


2.2 Landau Fermi liquid and its topological stability 


Let us start with the Fermi gas - the system of free fermions with the spectrum 
gp = p^llm — p, where ju >0. The energy functional for free fermions is ^{n(p)} = 

/ EpWp, which gives the solution shown in Fig. j^(/e/f). The solution of the 

class gp = 0 forms the Fermi surface with radius pp, where p^/2m = p. Outside of 
the Fermi surface the distribution function rip = const, with «p = 1 for p < p/r and 
lip = 0 for p > pp. This corresponds to the class of solutions with drip = 0. 

The structure with the Fermi surface is topologically protected. The topological 
stability can be illustrated using the d = 2 system, when the Fermi surface p = 
Pp is the Id circle. For that one has to consider the Green’s function at imaginary 
frequency: 


G(a),p) 



( 11 ) 
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Fig. 2 Illustration of the topological stability of the Fermi surface on an example case with dimen¬ 
sion d = 2, when the Fermi surface forms a closed loop. Green’s function has singularities on the 
line CO = 0, pj + Py = Pf in the three-dimensional space (a),p^,py). Stability of the Fermi surface 
is protected by the invariant GD which is represented by an integral over an arbitrary contour C 
around the Green’s function singularity. 


The Green’s function in Eq. (Ill has singularities at O) = 0 for p belonging to the 
Fermi surface. These points form a closed line in the three dimensional {co,px,Py)- 
space, see Fig. |^for d = 2. This line has a topological winding number; the phase 
<P of the Green’s function, G = |G|e‘*’ changes by 2n along an arbitrary contour C 
around the line. In other words, the Fermi surface represents the p-space analog of 
the vortex lines in superfluids and superconductors, where the phase of the order 
parameter changes by 27Z around the vortex. The 2n winding of the phase <P can¬ 
not change under small deformations of the parameters of the system, and thus is 
robust to the interactions between the particles, if we do not consider the supercon¬ 
ducting, magnetic or other phase transition, which drastically (non-perturbatively) 
reconstructs the energy spectrum. This topological stability is the reason why inter¬ 
acting Fermi liquids preserve the Fermi surface. 

For more complicated cases, when the Green’s function has spin, band and other 
indices, and for arbitrary dimension d the winding number N of the Fermi surface 
is expressed analytically in terms of the matrix Green’s function in the following 
form; 


N = tr j^^^.Gico,p)d,G-\co,p). (12) 

Here the integral is taken over an arbitrary contour C around the Green’s function 
singularity in the d+ \ momentum-frequency space. For the Green’s function in 
Eq.([TT) the topological invariant = 1. 


2.3 Khodel-Shaginyan flat band and its topology 

When the interaction between the particles is strong enough, so that it starts dom¬ 
inating over the fermionic statistics, a more classical behavior of the distribution 
function may emerge, in which the natural solution of the variational problem corre- 
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Fig. 3 Illustration of the topological structure of the Khodel-Shaginyan flat band in a = 2 system. 
According to Fig. ^{right) the Fermi surface al p = pp spreads into a flat band concentrated in 
the region pi < p < p 2 - Correspondingly the line of the Green’s function singularities in Fig.[^- 
a vortex line - is spreading to an analog of a domain wall terminating on a pair of ;r-vortices at 
p = Pi and p = P 2 - (H 

sponds to a zero value of the variational derivative, 5S' {n(p )}/drip = 0. Therefore, 
with an increasing interaction strength one may expect the topological quantum 
phase transition to the distribution in Fig. {right), where the solution gp = 0 is 
spread over a finite region in momentum space, i.e., forming a flat band. 

The topological structure of the Khodel-Shaginyan flat band is shown in Fig. 
The vortex line with 2n winding transforms to the domain wall, at which the Green’s 
function has a jump, G{(0 = +0,px,Py) — G{(0 = —0,px,Py) 7 ^ 0. This domain wall 
terminates on 7Z vortices. 


2.4 Flat band near the saddle point 

Here we consider the formation of the Khodel-Shaginyan flat band in the vicinity 
of a saddle point in the d — 2 quasiparticle spectrum using the phenomenological 
Landau theory |7] and compare it with the results of numerical simulations of the 
Hubbard model [81 . 

The simple Landau-type functional for the interacting Fermi liquid is 


+ T L /(P>P')«P«p ■ 

D n n' 


(13) 


P ^ p.p' 


We illustrate the flat band solution using an even simpler functional with contact 
interaction; 



(14) 


where U >0. This functional has always a flat band solution with 0 < Wp < 1; 



( 15 ) 








Flat bands as a route to high-temperature superconductivity in graphite 


7 




Fig. 4 Flat band emerging near a saddle point. {Left): from the simplified Landau-type theory in 
Eqs. jl4^ -l|16[l with [1 = 0. The flat band is concentrated in the black region. {Right): from the nu¬ 
merical solution of the Flubbard model |8( [D. Yudin, et al, Phys. Rev. Lett. 112, 070403 (2014)], 
showing the spectral function within the reciprocal space of an interacting triangular lattice. The 
lower left sextant corresponds to the noninteracting case U = 0. For large U the band flattening is 
clearly seen near the saddle points. 


«P = 2- 


U 


0 < Hp < 1 . 


( 16 ) 


In the vicinity of the saddle point the non-perturbed spectrum (i.e. at 1/ = 0) has 
the form — fi. For /i 0, there are two hyperbolic Fermi surfaces. They 

interconnect at the Lifshitz transition, which takes place at ji =0. When the inter¬ 
action U is switched on, the flat bands emerge. Figuredemonstrates the flat 
band at jj. =0. The same shape of the region with the flat band has been obtained 
from the numerical simulations of the Hubbard model lO, see Fig.^(riglit). 


3 Flat bands in topological semimetals 

In the previous section we discuss the fermion condensate - the flat band, which 
may emerge due to interactions in the vicinity of a singularity in the non-interacting 
spectrum. There are also other ways to generate flat bands or approximate flat bands. 
Historically the flat bands first appeared as Landau levels of charged particles in 
a magnetic held ||9l. Here we discuss the flat bands that have purely topological 
origin. They may exist without a magnetic held, and they are not very sensitive to 
interactions. The flat band may emerge as the surface or interface state in topological 
semimetals Col [mill, which we discuss in this section. The flat band may also 
appear at the strained interfaces with misfit dislocations, which play the role of 
effective magnetic held. Such flat bands are discussed in Sec.j^ 

In this section we characterize semimetals that have an internal spin-like struc¬ 
ture. All the examples discussed here can be characterized via the topological in¬ 
variant of the form 
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Ni=tr£^.-[rH-\p)d,H{p)l (17) 

where H is the Hamiltonian in the momentum space; F is a matrix which commutes 
or anti-commutes with the Hamiltonian, such as the third Pauli matrix ff 3 acting on 
the spin-like degree of freedom; and C is a contour in momentum space, specified 
separately for each semimetal. 

In semimetals, the flat bands are realized as a particular consequence of the bulk¬ 
boundary correspondence of topological media (see lfT4l [TSl [Thl or Sec. 22.1 in 
ini). If for example a 3D bulk system contains Weyl points, then an interface with 
a topologically trivial material, or with a material having a different value of topo¬ 
logical invariants contains a line of zeroes - the Fermi arc m. The termination 
points of the Fermi arc are given by projections of the Weyl points to the interface. 
In the same manner the Dirac lines in 3D bulk or Dirac points in 2D bulk give rise 
to the to nodes of higher dimension at the interface - the flat band ifTOlfTTlfT^ . The 
boundaries of the flat band are determined by the projection of the Dirac line or 
Dirac points to the interface. 


3.1 Topological nodes: Dirac lines and Dirac points 

Let us consider an example semimetal characterized by an even number of 2- 
dimensional Dirac points, such as that found in graphene around the two val¬ 
leys. Close to the Dirac point the Hamiltonian can be written as a 2x2 matrix in 
(pseudo)spin space. 


H = vf 


0 

Px + ipy 


Px - ipy\ 

0 J 


= vfP • a, 


(18) 


where p = {px,Py) and a — ((7^,0).), where Gj are Pauli matrices. The eigenvalues 
of H satisfy having a node e(p = 0) in a single point in momentum 

space. To illustrate the topological protection of this node, let us add a perturbation 
of the form V (p) = v(p) • (J, where v = {yx,Vy) does not break the (pseudo)spin sym¬ 
metry. As a result, the dispersion becomes = {vpPx + Vxf' + {vppy + Vy)^. which 
again has a single node at {px,Py) = —{vx,Vy)/vF- The only effect of the potential 
is thus to shift the node, but not annihilate it. This property can be expressed via the 
presence of the topological charge of the form ([T^, where F = G^ and the contour C 
goes around the Dirac point in the 2D momentum space (see Fig.[^. For the Hamil¬ 
tonian in Eq. we get Ai = 1. In graphene, there are two Dirac points: the first 
one has the form in Eq. ( fTS) !, and the second is otherwise the same but py i— —Py- In 
that case the second Dirac point has Ni = — 1. These topological charges stay invari¬ 
ant to perturbations of the form V (p), as long as we shift the contour of integration 
along with the shift of the Dirac point, and as long as the two nodes do not merge 
due to such a perturbation. 
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Fig. 5 2D Dirac/Fermi point in momentum space and the line of integration for the topological 
invariant Ni. 

In 3D materials, the topological charge in Eq. ([T7J characterizes lines of nodes 
- the Dirac lines. The Dirac lines are readily obtained in superfluids and supercon¬ 
ductors, where the symmetry operator F contains the particle-hole symmetry. In 
particular, the Dirac line exists in the polar phase of superfluid ^He ifT^ and may 
appear in superconductors without inversion symmetry m- The topologically sta¬ 
ble Dirac lines give rise to the topologically protected surface flat band. According 
to the bulk-surface cotTespondence, the boundary of the flat band is determined by 
the projection of the nodal line on to surfaces. In nonsuperconducting materials the 
corresponding symmetry which enters Eq. ([T^ can be only approximate, being vio¬ 
lated by spin-orbit interaction, or by the higher order hopping elements. This leads to 
formation of approximately flat surface bands as it happens for example in graphite 
ll20l . graphene networks ED and possibly in some other materials ll22l . 



Fig. 6 Energy spectrum in the conduction and valence bands of graphene. Flere a is the lattice 
constant, and yi denotes the nearest-neighbour hopping parameter in the tight-binding lattice. The 
underlying contour plot shows the positions of the Dirac points. Only two of these points are non¬ 
equivalent, the others are connected via reciprocal lattice vectors. 









10 


T.T. Heikkila and G.E. Volovik 


We first consider the case of graphene, whose energy spectrum within the valence 
and conduction bands is plotted in Fig. This spectrum follows for example from 
the nearest-neighbour tight-binding model on a honeycomb lattice, see Enin. 
Around the specihc points, marked K and K' in the hgure, the low-energy Hamilto¬ 
nian is of the form of Eq. ( |TF| . These points are described by the topological charge 
Ni, such that Ni (K) = -|-1 and Ni {K') = —1. 

Figure]^ shows the formation of flat bands at the edges of graphene ifTOll . There, 
we show the locations of the K and K' points in the 2d momentum plane, and 
consider the presence of an edge placed in the .r-direction. Now, the edge marks 
a boundary between graphene, which is a nodal semimetal, and vacuum, which is 
a trivial insulator. By the bulk-boundary correspondence we may hence expect flat 
band states at the edge. Since we maintain translational invariance along the x di¬ 
rection, Pjc remains a good quantum number and it also parametrizes the edge states 
and their dispersion Eeipx)- According to the bulk-boundary correspondence, the 
projections of the Dirac points to the boundary determine the end points of the flat 
band. From Fig. [7] it is clear that the termination points of the flat band must be lo¬ 
cated at px values corresponding to the pj:-component of the K and K' points. When 
Pjf crosses these points, the normal to the interface run across either K or K' points 
and the topological invariant in Eq. ( [T7| ) along the normal changes. However, from 
the figure alone one cannot say which of the regions between the K and K' points 
contain flat bands and which not. The solution is to construct the system by repeat¬ 
ing a set of infinite chains, and construct the topological invariant, in this case called 
the Zak phase ll^ for these chains. As shown for example in l^fToll . the details 
of this procedure depend on the microscopic form of the edge, which is either of the 
“bearded” or the “zigzag” type. Figure [^illustrates the resulting positions of the flat 
band dispersion in momentum space. 

If we were to place the interface in the y direction (in which case we would obtain 
the “armchair” edge), there would be no flat band due to the fact that the normal to 
the interface would run across both K and K' points. In this case the change of the 
topological charge across this normal line would be Ni (K) +Ni {K') = 0. 

In the following sections we explain the topological invariants in multilayer 
graphene structures and discuss the consequent (approximate) flat band surface 
states. 


3.2 From graphene to graphite 

In graphite, the coupling between the graphene layers is much weaker than that 
within the atoms producing the hexagonal lattice in individual layers. Because of 
this, we may assume that the individual layers remain to be described by the two- 
dimensional low-energy momentum-space Hamiltonian of the form in Eq. ( [TSl l. In 
the (most stable) AB stacking, the pair of layers is arranged so that the graphene 
hexagons are rotated with respect to each other by 30 °. As a result, one pair of 
atoms (say, A on the bottom layer and B on the top layer) reside on top of each 
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Fig. 7 Formation of a flat band in the zigzag (lower part of the figure) or bearded (upper part) 
edges of graphene. 


Other, whereas the other atom of the unit cell (B on the bottom and A on the top) are 
on the bottom/top of the other hexagon. Therefore, the interlayer coupling between 
the first pair of atoms is much stronger than the coupling between the second pair. 
Taking into account only this strongest coupling then produces the Hamiltonian of 
the form 


Hk 


VfP • a 


-rio’i 

VfP • a 


(19) 


around one of the graphene valleys (/f-point) for this pair of layers. Here 71 quan¬ 
tifies the interlayer hopping. Around the other valley {K' point) the Hamiltonian 
is Hi^t{px,Py) = Hxipxj—Py)- It is straightforward to show that have four 

branches of eigenvectors, two gapped ones (with e(|p| = 0 ) = ± 7 i) and two with 
a quadratic dispersion around a Fermi node, = \Yi?'{pIPfbY, where ppB = 
| 7 i l/vp. The corresponding N\ = 2 (—2) for Hk {Hk')- 

Beyond the bilayer, the graphene layers can be stacked in two qualitatively differ¬ 
ent ways by respecting the AB stacking for each pair of layers. In Bernal stacking, 
the line of strongest interlayer coupling is straight, i.e., connects the same atoms in 
each layer whereas for rhombohedral (ABC) stacking, it follows an armchair-type 
pattern, i.e., the strongly coupled atoms differ between neighbouring pairs of layers. 
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The corresponding elements of the momentum space Hamiltonians coupling layers 
n and m for 2D momenta around the K point are 

vfP • a5mn - 7i5m.«+i {oy[(- 1)'" - 1] + + 1]} /2 + h.c. 

( 20 ) 

fFP'Cr^mn 7l ^H?,n+1 T h.C. (21) 

In the terms in square brackets take care of indexing the even and odd layers 

separately. Note that since the choice of the A /B indices for the graphene sublattice 
atoms is arbitrary, we could have also chosen to write the coupling in //rhombohedrai 
in terms of (J^ instead of (J;. Formally, these two choices correspond to different 
signs of the bulk topological invariants (see below), but they have a meaning only 
in the presence of (Bernal) stacking faults that change the coupling around some 
particular interface. 


H. 


Bernal 


ti. 


rhombohedral 


3.3 Bulk Fermi line in Bernal graphite 


We first derive the bulk dispersion of Bernal graphite by including only the strongest 
interlayer hopping term y\. We make the bulk Ansatz 




(So) n odd 

«£. j3e) , n even 


( 22 ) 


for the wave function on the n-th layer. Here is the momentum in the direction 
perpendicular to the layers and a is the distance between the layers. This is an eigen¬ 
function with energy e provided that the coefficients satisfy 


/ 0 

vfP- 

0 

- 27 icos(p,a/;j)\ 

fa,\ 


/aA 

VFP+ 

0 

0 

0 


Po 

- £ 

Po 

0 

0 

0 

vfP- 


Ue 


ae 

^- 27 ^ cos{p^a/h) 

0 

VFP+ 

0 

/ 

\Pe) 




Hb 


where p±= ipy = p±e“<’, p± > 0. This has (four) zero-energy solutions at 
Px = Py = 0, regardless of the value of p^. Within this approximation, Bernal 
graphite has therefore two Dirac lines, running through the K and K' points of the 
2D graphene band structure, or between the H points in the 3D graphite band struc¬ 
ture, see Fig. 3.3 Including higher-order hopping terms in the Hamiltonian expands 
this line into electron and hole pockets ll26l . However, let us first analyze the topol¬ 
ogy of Hg. The topological charge in this case is of the form of Eq. ( [T7| ), where 
the contour C runs around the Dirac lines as indicated in Fig. 8, and ffj should be 
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replaced by (ffj.® 1). This produces A^i = ±2, where the factor 2 takes care of the 
additional layer degree of freedom in Eq. (|23|). 



Fig. 8 Fermi lines in Bernal graphite within the nearest-neighbour interlayer approximation. This 
picture expands the 2D momentum space (pxjPy plane) of graphene to the 3D space where runs 
between —Tin!a to fitr/o (a is the interlayer spacing). 


This is thus nothing but the generalization of the graphene topological charge to 
graphite and it yields = ±2 for the momenta along the lines H — K — H/H' — 
K' — H where the point is shifted from the K^''^ point in the direction 
by hn/a. Based on the bulk-boundary correspondence we may therefore expect to 
have surface states at the lateral boundaries of Bernal graphite, as extensions of the 
flat band states in zigzag graphene, but now the flat band extends throughout the 
first Brillouin zone in the direction. 

Let us now consider the effect of higher-order hopping elements. First they split 
the Dirac line with multiple charge N\ =2 into elementary Dirac lines. The natural 
splitting would be into two lines with charges = 1. However, the situation is 
more interesting, see Fig. I^and Refs. The N\ =2 line splits into 4 Dirac 

lines: three lines with trigonal arrangement have A^i = 1, while the central line has 

= — 1. The total topological invariant remains the same, A^i = l-|-1-|-1 — 1=2. 
In addition, the higher-order hoppings break the symmetry F EOl underlying the 
topological protection. Without the topological protection the four Dirac lines are 
expanded into four pairs of electron and hole pockets, which touch each other at 
four points. As a result the surface states of Bernal graphite cease to be flat on the 
energy scale related to those higher-order hoppings. 


3.4 Spiral Fermi line in rhombohedral graphite 

To obtain the bulk dispersion of rhombohedral graphite, it is enough to make the 
Ansatz 
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electrons 



Fig. 9 Graphite Fermi surface from Ref. (26l . It consists of the electron and hole pockets connected 
at four points. Each point is described by topological invariant in Eq. CZ): the point in the center 
has Ni = — I and the other three points have Vi = 1. Such structure originates from the nodal 
line with Vi = 2 in the Hamiltonian l |23| l, when the higher order hopping elements are taken into 
account. 



Wn = (24) 

This is an eigenfunction of //rhombohedrai ^ satisfy 

f 0 VFp^-Yiexp{ipM/h)\fa\_ fa\ 

\vFp+-f^exp{-ipM/n) 0 jyp) \p) 

" -V-^ 

Hrhg 
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for some energy £. The zero-energy solutions are found along a spiral in the mo¬ 
mentum space. 



(26) 


In this case the topological charge equals the spiral helicity, and can be defined by 
Eq. ([T7]) where the integration contour runs over the 1st Brillouin zone (—nja to 
%ld) in the p^-direction. For running close to the line from the H point via K 
back to H, = 1 whenever the transverse momentum p^i /p, is inside the spiral, 
i.e., px < |7i l/vf. For larger momenta away from the H — K — H line, = 0. On 
the other hand, for the path H' —K' — H' the spiral is defined by Eq. ( [26l ), where 
Py O —Py, inverting the helicity and the sign of A^i to -1. 

As a result of the existence of the non-trivial topological charge, the surface 
layers of rhombohedral graphite contain flat bands, i.e., e(px < |7i|/vx) = 0. This 
can be understood as follows. For each px with px ^ |7i|/''x, the Hamiltonian 
H^yiPz) describes a ID insulator. For px < |7i|/vx, this insulator is topological, 
since it has the topological charge Ati(px) = 1. According to the bulk-boundary 
correspondence, each topological insulator has an edge state with zero energy. These 
states form a flat band in the region px < |7i \/vf- 


4 Flat band at strained interfaces 

Recently another possible source of the topological flat band has been discussed 
in two materials: highly oriented pyrolytic Bernal graphite (HOPG) 12^ and het¬ 
erostructures SnTe/PbTe, PbTe/PbS, PbTe/PbSe, and PbTe/YbS consisting of a topo¬ 
logical crystalline insulator and a trivial insulator 1^ . In both cases the flat band 
comes from a misfit dislocation array, which is spontaneously formed at the inter¬ 
face between two crystals due to the lattice mismatch. In Ref. Il28ll the lattice of 
screw dislocations has been considered, which emerges at the interface between 
two domains of HOPG with different orientations of crystal axes. In Ref. Il2^ the 
misfit dislocation array is formed at the interface between topological and trivial 
insulators. 

The above two systems exhibit similar phenomenon. In both cases superconduc¬ 
tivity related with the interfaces has been found Il30l [311 . The reported transition 
temperature essentially exceeds the typical transition temperature expected for the 
bulk materials. A possible origin of this phenomenon is the flat band at the inter¬ 
faces, where the transition temperature could be proportional to the coupling con¬ 
stant and the area of the flat band. 

The topological origin of the flat bands in these systems can be understood either 
in terms of the overlapping of the ID flat bands formed within the dislocations 
or using the following consideration. In case of the heterostructures, on one side 
of the interface the insulator is topological, and thus the interface contains Dirac 
fermions. The strain at the interface (with its dislocations) acts on Dirac fermions 
as the effective magnetic held. Such emergent held is now extensively discussed for 
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Fig. 11 Misfit dislocation grid at the interface from IN. Ya. Fogel, Phys. Rev. B 66, 174513 
(2002)] (30). 


Strained graphene, see e.g. the most recent paper Ref. If3^ and references therein. 
The effective magnetic field B produces the required flat band, since the first Landau 
level for massless Dirac fermions has zero energy. When the period decreases, the 
field B increases, giving rise to enhanced density of states, which is proportional to 
B. This is equivalent to increase of the area of flat band in the scenario discussed in 
Sec.[^ Such increase of transition temperature is discussed in Ref. 1291 . 

In particular. Ref. l29l considered the case where the massless Dirac Hamilto¬ 
nian of the topological insulator experience a periodic field. They arrived at the 
Hamiltonian of the form 


Htf = -idxOx + -3l W]cTv, (27) 

written in terms of (scaled) momentum ky along the dislocation, and a scaled strain- 
induced gauge field A(x) satisfying A (x-ft/) = A(x) and Jq A{x)dx = 0. In the case 
of Ref. l29l . A(x) = j3 cos{2nx/d), but the approach works for a more general pe¬ 
riodic vector potential as well. Due to the position dependent vector potential, there 
is no translation symmetry in the x direction and the momentum in the x-direction 
is not a good quantum number. However, due to the periodicity of the potential we 
can use Bloch’s theorem and define the pseudomomentum kx- This allows for calcu¬ 
lating the spectrum of Eq. plotted in Fig. and exhibiting (an approximate) 
flat band for ky <15/2 and for all values of kx G [—n/d,7t/d]. 

The emergence of the approximate flat band in Eq. ( |27l l can be understood by first 
considering the case ky = 0. In that case the Hamiltonian has two (unnormalized) 
zero-energy solutions. 
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ky 


Fig. 12 Spectrum of Eq. with P = 30. Inset shows the (approximate) linear dispersion with 
low values of ky, and with a speed as in Eq. j^. 


¥+ = exp[^ A{x')dx'] 

V/- = A{x )dx]. 


Let us include the term H\ = kyOy as a perturbation. The second order secular equa¬ 
tion produces the 2x2 Hamiltonian constructed from the matrix elements of the 
Hamiltonian H \; 





with the spectrum E = -Ecky. The slope c of the spectrum is obtained from the matrix 
elements of the normalized eigenfunctions 


/“^£/xexp[2 /q A(x')c/ji:'] e^IoMx'jdx' 


For A{x) = j3cos(x) we thus get c = l//o(2j3), where /o(x) is the zeroth Bessell 
function of the first kind. When c ^ 1, requiring j3 7^ 1, the result is an approximate 
flat band. 

Alternatively, we may consider the case of a periodic line dislocation. There, the 
vector potential is proportional to the strain field, A(x) = aiU;yx{x) + a 2 Uyy(x), where 
a\ and a 2 are coupling constants and 
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^ ^ _ bz 3x^+z^ 

“ 2k{1-v) (x2+z2)2 

, ^ _ bzv 1 

7r(l-v)x2 + z2’ 


where b is the size of the Burgers vector of the dislocation, z is the distance from the 
dislocation plane, and v is the Poisson ratio. We also assume that such dislocations 
repeat after each d. In this case we may estimate the slope c for ^ z by 



[ r At 


r b{ai+va2)'\ 

c ~ exp 

—2 / Ahc) 
Jo 

= exp 

1-v 


The flat band dispersion where c <C 1 thus requires either v « 1 or coupling con¬ 
stants a, ^ 1 jb. Note that this result show that the estimates done in Ref. Il29ll are 
overoptimistic, as there ab~\. 


5 Superconductivity in graphene or graphite? 

As discussed in Sec. bulk neutral graphite and graphene are typically considered 
to have a very low density of states. This is why the occurrence of superconductivity 
in these bulk systems would require a strong doping, shifting the Fermi energy to a 
finite value and thereby increasing the density of states. Let us consider the example 
case of a 2D Fermi line with the superconducting coupling g, but with a shifted 
Fermi level by a value ji. In this case Eq. Q yields for ji —0 ll33l 

A = (fle„v-l/a)/2>0, (28) 

where a = g/{n'h^v\). This is non-zero only for a large coupling strength a >1/ 

For weak coupling |fl| < 1 superconductivity is enhanced by doping, because for 
non-zero jx the solution is A = 2|/r|exp(— ~ !)■ 

This strategy has been followed by some graphite experimentalists. For example 
Ref. 041 demonstrates bulk superconductivity of two graphite intercalation com¬ 
pounds CeYb and CeCa with critical temperatures of 6.5 and 11.5 K, respectively. 
There one of the effects from the intercalant layers is the charge transfer, i.e., doping 
the graphite. 

The situation is opposite in flat band systems; doping does not aid supercon¬ 
ductivity, but rather decreases the critical temperature. Fith a non-zero chemical 
potential /r, the gap becomes 


A = 



( 29 ) 


where Afb is the gap obtained for /r = 0, as in Eq. 0 - 
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There are experimental indications for the presence of superconductivity at 
graphite interfaces El. In this case, the interfaces form between differently ori¬ 
ented regions of Bernal graphite, so that the graphite c-axis (p^ in our work here) 
runs perpendicular to these surfaces. These observations hence cannot be explained 
by the flat bands at the lateral interfaces (along px or py) as one would expect in 
Bernal graphite. There are also no direct indications of the presence of rhombo- 
hedral stacking in these systems, but the resolution of the imaging does not allow 
ruling out some rhombohedral regions close to the interfaces exhibiting supercon¬ 
ducting properties. Note that rhombohedral stacking has been found in highly ori¬ 
ented pyrolytic graphite, see for example ll35l[^ . The remaining scenario is based 
on the formation of an array of dislocations (Sec.|^ at these surfaces as the surface 
layers try to adapt to the neighbouring graphite planes. 

The discussion in Sec. [^concerns only the relation of the mean field order pa¬ 
rameter with normal-state electronic spectrum. The details of the flat band supercon¬ 
ducting state, such as the quasiparticle spectrum, supercurrent and collective modes, 
depend on the way the flat band is formed. The case of rhombohedral graphite has 
been worked out in Refs. IIJTIIMII^ and summarized in ll40l . The density of states 
in the superconducting state was computed in m- In that case the superconducting 
spectrum is not flat, but rather has an inverse parabolic shape and exhibits a minigap 
whose value is inversely proportional to the distance between adjacent surfaces. The 
supercurrent is characterized by a large critical current that is linearly proportional 
to 4. However, in this model the supercurrent does not only flow along the surfaces, 
but also between them. On the other hand it is clear that not all types of flat bands 
can describe a state with a non-vanishing critical current. For example, an entirely 
isolated region in momentum space having a flat band describes a set of completely 
localized electrons, which cannot be made to carry supercurrent. 

In this chapter we have concentrated on illustrating the results for mean-held 
models of superconductivity on flat bands. As discussed in HOt . fluctuation effects 
are much more important in flat band superconductors than in ordinary BCS super¬ 
conductors. The nature of the fluctuation response depends on the mechanism of the 
flat band formation, and the ensuing superconducting spectrum. The details of such 
fluctuation effects we leave for further work. 

As we have discussed in this chapter, there are many ways of realizing (approxi¬ 
mate) flat bands. Some of these realizations may be relevant in interfaces of graphite, 
where indications of superconductivity with a very high critical temperature have 
been observed. It remains to be demonstrated whether these properties can be ex¬ 
plained within the flat band scenario of superconductivity. 

We acknowledge Pablo Esquinazi, Ville Kauppila and Timo Hyart for discus¬ 
sions. This work was supported by the Academy of Finland through its Center of 
Excellence program, and by the European Research Council (Grant No. 240362- 
Heattronics). 
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